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Abstract 

The NASA Solar TErrestrial RElations Observatory (STEREO) mission opencd| 
up the possibihty to forecast the arrival times, speeds and directions of solar 
transients from outside the Sun-Earth line. In particular, we are interested 
in predicting potentially geo-effectivc Interplanetary Coronal Mass Ejections 
(ICMEs) from observations of density structures at large observation angles from 
the Sun (with the STEREO Heliospheric Imager instrument). We contribute to 
this endeavor by deriving analytical formulas concerning a geometric correction 
for the ICME speed and arrival time for the technique introduced by Davies et 
al. (2012, Astrophys. J., in press) called Self-Similar Expansion Fitting (SSEF). 
This model assumes that a circle propagates outward, along a plane specified 
by a position angle (e.g. the ecliptic), with constant angular half width (A). 
This is an extension to earlier, more simple models: Fixed-$ Fitting (A = 0°) 
and Harmonic Mean Fitting (A ~ 90°). This approach has the advantage that 
it is possible to assess clearly, in contrast to previous models, if a particular 
location in the heliosphere, such as a planet or spacecraft, might be expected 
to be hit by the ICME front. Our correction formulas are especially significant 
for glancing hits, where small differences in the direction greatly influence the 
expected speeds (up to 100 — 200 km s~-^) and arrival times (up to two days later 
than the apex). For very wide ICMEs (2 A > 120°), the geometric correction 
becomes very similar to the one derived by Mostl et al. (2011, Astrophys. J. 741, 
id. 34) for the Harmonic Mean model. These analytic expressions can also be 
used for empirical or analytical models to predict the 1 AU arrival time of an 
ICME by correcting for effects of hits by the flank rather than the apex, if the 
width and direction of the ICME in a plane are known and a circular geometry 
of the ICME front is assumed. 
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1. Introduction 

Predicting the effects of solar eruptions at Earth is a scientific art in the making. 
Many different approaches have been developed world-wide, using analytical, 
numerical or empirical techniques, to better understand the physics behind 
Interplanetary Coronal Mass Ejections (ICMEs) which should ultimately lead 
to space weather predictions on a routine basis with reasonable accuracy. The 
NASA Solar TErrestrial RElations Observatory (STEREO, Kaiser et al. 2008), 
launched in 2006 and consisting of two observatories leading and lagging the 
orbit of Earth by 22.5° per year, plays a major role in these efforts. STEREO 
has been designed to enhance our ability to forecast ICME arrival times at Earth 
through more precise determination of ICME speeds and directions with both 
stereoscopic and single-spacecraft techniques, based on observations taken from 
outside the Sun-Earth line to get rid of projection effects (Gopalswamy et ai, 
2001; Schwenn et ai, 2005) and filling the imaging gap between the Sun and 
1 AU with the Heliospheric Imagers (HI) instruments (Eyles et al., 2009). The 
latter has also been achieved from within Earth orbit by the Solar Mass Ejection 
Imager (SMEI) instrument (Eyles et ai, 2003). 

An idea originally introduced by Shcclcy et al. (1999) was used to predict 
directions and speeds of ICMEs (Davis et ai, 2009) and co-rotating interaction 
regions (Rouillard et al. 2008) from HI images. This method makes use of a 
simple geometrical fact, which can be illustrated as follows: if someone throws 
a ball with constant velocity along a constant direction, an observer measuring 
the angle between himself, the thrower and the ball as a function of time will 
notice deceptive decelerations or accelerations in the time-angle profile. The 
form of this function depends on the speed and direction of the ball with respect 
to the observer. Applied to ICMEs in any plane given by a position angle to 
solar north (e.g. the ecliptic plane), these profiles, known analytically for either 
point-like (Fixcd-'J, FP, Sheeley et al. 1999; Kahler and Webb, 2007) or circular 
fronts (Harmonic Mean, HM. Lugaz, Vourlidas, and Roussev, 2009; Howard and 
Tappin, 2009), can be inversely fitted through a minimization process to observed 
profiles of the elongation angle of the ICME front from the Sun (e.g. Rouillard 
et al. 2009; Davis et al. 2009; Howard and Tappin, 2009; Tappin and Howard, 
2009; Davies et al. 2009; Savani et al. 2009, 2010; Mostl et al. 2009, 2010, 2011; 
Lugaz, 2010). We denote the techniques as EPF and HMF. The EPF technique 
has been successfully applied in real time to an ICME in April 2010 (Davis 
et ai, 2011), and to some other recent real time events, and is a candidate 
for a technique to be used routinely on a future space weather mission at 
the L5 point in the Sun-Earth system (Gopalswamy et ai, 2011) or on Solar 
Orbiter. In comparison to a numerical simulation, some biases have been found 
for FPF/HMF, mainly for cases where the ICME is either close to the limb or 
heading towards the observer (Lugaz, Roussev, and Gombosi, 2011). Similarly 
to the above mentioned techniques, stereoscopic versions exist for FP (Liu et al. 
2010a) and HM (Lugaz et ai, 2010), relaxing the assumptions of constant speed 
and direction, though these techniques have not yet revealed strong deviations 
from radial propagation far from the Sun. 
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An improvement to the single-spacecraft techniques has been introduced by 
(Davies et at 2012, further cahed DA12), called self-similar expansion or SSE 
to convert elongation to distance. This is based on model 2 by Lugaz et al. 
(2010); and DA12 made it suitable as a prediction tool through inverse fitting 
(SSEF). This model has a great appeal by solving some problems (see discussion 
by Mostl et al. 2011, further: MOll) which are intrinsic to the point-like or very 
wide front assumptions of FP and HM, respectively, especially when thinking 
about when and with which speed the ICME front will hit a planet, resulting 
in possible magnetic storms, or a spacecraft which makes in situ observations of 
the solar wind plasma and magnetic field parameters. 

The aim of this short report is to clarify the calculation of arrival times and 
speeds at a particular position in the heliosphcre with the self-similar expan- 
sion model (DA12) or its stereoscopic version, model 2 in Lugaz et al. (2010). 
The assumption behind it is that the ICME front can be approximated by a 
circle with constant angular width (see Figure 1). Our term "front" leaves open 
whether one wants to investigate the speed and arrival time of the interplanetary 
shock wave or the front boundary of the magnetic flux rope or ejecta driving 
the shock. Given the position of a spacecraft with respect to the direction of the 
ICME apex (the point along the front with greatest heliospheric distance from 
the Sun), this model will lead to later arrival times and lower speeds, or even no 
hit at all. In this way, we extend the work of MOll who derived a similar arrival 
time correction for the HM model, which is a limiting case of SSE for large ICME 
width. The FP model can also easily be understood as the limiting case of SSE for 
negligible width. Clearly, in the real solar wind environment it is questionable if 
these idealized conditions are even roughly met, and distortions of the fronts will 
be likely. Also note that a spherical expanding front centered on the Sun (called 
"Point-P" by various authors) has the same heliocentric distance everywhere, 
and corrections to speed and arrival time are not needed. However, to derive 
geometrically consistent arrival times and speeds with SSE or its triangulation 
version it is necessary to use the formulas derived in this paper. 

The analytical considerations we present in this paper are not restricted to 
techniques used on HI instruments, but they provide a simple framework to 
assess the effects of an ICME hitting a target with its apex or flank. If the 
direction of an ICME is known from stereoscopic observations close to the 
Sun from various techniques (Thernisien, Vourlidas, and Howard, 2009; Tappin 
and Howard, 2009; Temmer, Preiss, and Veronig, 2009; de Koning, Pizzo, and 
Biesecker, 2009; Mierla et al, 2010; Byrne et al, 2010; Rodriguez et al, 2011), 
an educated guess can be made regarding the difference in heliospheric longitude 
between the ICME apex and the spacecraft or planet. A model for ICME prop- 
agation can then be used to extrapolate the arrival time and speed of the ICME 
apex to 1 AU (e.g. Gopalswamy et al 2001; Vrsnak and Gopalswamy, 2002; 
Schwenn et al 2005; Siscoe and Schwenn, 2006; Vrsnak et al. 2010; Maloney 
and Gallagher, 2010; Vrsnak et al 2012, this issue), which can subsequently be 
corrected for the longitude difference using the formulas presented here, given 
that the width of the ICME is known, or simply assumed. 

The average width of a CME in coronagraph observations is known to be 
about 50° (Yashiro et al, 2004) with considerable scatter around this value. 
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Direct comparison of STEREO/HI movies to a rendered 3D density model of 
a loop-like ICME revealed a true (de-projected) width of 70° (Wood, Howard, 
and Socker, 2010) for one event. However, we are interested here in the angular 
extent of the ICME front along the line-of-sight (LOS) in a given plane (such as 
the ecliptic). Stereoscopic modeling assuming such a geometry by Lugaz et al. 
(2010) revealed half (fuh) widths in the ecliptic from 25 - 45° (50 - 90°), and 
observations of interplanetary shocks and their solar sources point to extents of 
up to 100° (Richardson and Cane, 1993). One should also be aware that the 
ICME angular extent along the LOS should depend upon its orientation with 
respect to the ecliptic plane (e.g. Liu et al. 2010b; Kilpua et al. 2011; Mostl et al. 
2011). 

We start by showing the derivation of the formulas in Section 2 and discuss 
them in comparison to the former FP and HM models. We then proceed to plot 
the effects on arrival time and speed predictions in the next two sections, discuss 
their significance and conclude with comments on our assumptions and future 
work. 



2. Derivation of arrival time and speed for SSE 

Figure 1 shows the geometry introduced by Lugaz et al. (2010) and elaborated 
by DA12. It is assumed that an observer looks along the tangent to a circular 
front propagating radially away from the Sun with constant angular width. The 
elongation-to-distance conversion for SSE can be derived as: 

doSin(e(<))(H-sinA) 

^SSE = . . ..s , ,^ , — ■■ — 7- (,lj 
sm(e(t) + <p) + sm A 

with do the radial distance of the observing spacecraft (called O) from the Sun, 
e(t) the measured elongation angle as function of time, A the angular half-width 
of the ICME, angle 4> the propagation direction as measured from the observer 
and RssE the heliocentric distance of the ICME apex. For the large and small 
triangles in Figure 1 we can write: 

sme{t) sin(7r — e — (/)) 



RssE + a do 

sin(7r - e - 0) = . (3) 

r + a 

We further need a relation between Rsse and r, the radius of the circle. Using 
RssE = r + c and r = csin A leads to: 

r{RssE, A = , 4) 

(1 + sm A) 

and can be understood as a radius of curvature, which is, however, coupled to 
the width A. It also is a function of the apex position, and thus time, because the 
circle expands as the front propagates outwards. Using Equation 4 in Equation 3 
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Figure 1. The geometry of the SSE model. An observer looks along the tangent to a circular 
front at an elongation e(t) to Sun center. The circle expands with constant angular half-width 
A along a constant ICME propagation direction (p. Note that the full width of the ICME is 2A. 



and solving for a, allows one to put a into Equation 2, and after some rearrang- 
ing, Equation 1 is obtained. For completeness, the circle's central heliocentric 
distance is: 

''"(1+sinA)- 

By comparison, the HM model assumes a circle which is attached to the Sun at 
all times (Lugaz, Vourlidas, and Roussev, 2009; Howard and Tappin, 2009), and 
this simply corresponds to Equation 1 for the case of A = 90°: 

_ 2do sin eft) 
sm[e(t) + (pj + I 

If the width is negligible (A = 0°), the FP conversion for elongation-to-distance 
follows from Equation 1 (Sheeley et al. 1999): 

dp sin eft) 
" sin(e(t) + 0) 

We now ask ourselves: how far does the apex of the ICME front need to travel 
away from the Sun until it hits the in situ spacecraft or planet 17 Clearly, for 
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a circular front this will be always greater than Ps heliocentric distance, di. As 
seen from Figure 2, the ICME apex is separated at an angle A to /, so we make 
the following ansatz: 



sin A sin 



a 



r c 
sin A sin(7r — a — A) 

r di 
Using r ~ c sin A in Equation 8 results in: 



(8) 
(9) 



/smA\ , , 

a = arcsm . (10) 

\ sm A J 

Substituting Equations 10 and 4 into Equation 9 leads to: 

,^ , ^ , _ rf, sinA(l + sinA) 

KissE{i^, A, djj - — —-, — ■ , . >s , . ^ . — r (.11) 

sm(arcsm(sin A / sm A j + A j sm A 

Using trigonometric identities, this can be simplified to: 

d, (l + sin A) 

i?.ss£(A, A, dO = '\ > -. (12) 

cos A + V sin A — sin A 

We designated this special distance as RissE, because it is the heliocentric 
distance of the ICME apex at the arrival time ta when the spacecraft is hit by 
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Sun 

Figure 3. Illustration of five SSE fronts vfith different angular half widths A. The in situ 
observing spacecraft is positioned at an angle A to the ICME apex. 



the front, as illustrated in Figure 2. Again, for HM with A = 90°, Equation 12 
indeed reduces to: 

RiHM ~ — —r, (13) 
cos A 

similar to the one used by MOll. For FP, one may assume that a small spherical 
front is centered around the assumed point, so simply Ripp = di. 

2.1. Arrival time 

If we further assume that the ICME front travels with constant speed Vsse 
for all times when Rsssit) > di, we can express the arrival time correction 
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tcSSE > as: 



tcSSE = T7 ■ (14) 

VSSE 



This formula states for how long the ICME apex has to travel, after it passed di , 
until / is hit by the front. Here, Vsse can also be seen as an asymptotic speed 
given by an ICME propagation model. 

For the SSE fitting technique (DA12), which results in constants for A, 0, 
Vsse, and the launch time tossE, the arrival time at the spacecraft / can be 
expressed as: 

, , Ri^ssiA^d^) 

taSSE = toSSE H 77 • (J-Oj 

VSSE 

It needs to be understood that the launch time tgssE is not determined from 
chromospheric or low coronal observations, but that it is rather a back-projection 
to the center of the Sun (where e(t) = 0) to obtain a time of reference to be 
able to calculate taSSE (see also MOll). It is assumed here that the speed Vsse 
is a constant for the range of elongation angles for which the time-elongation 
track, which is fitted with SSEF, is extracted from HI observations. Any earlier 
accelerations or decelerations of the CME closer to the Sun do not affect the 
calculation of tassE- However, the launch time tosS-B will consequently be only 
a rough estimate of the real launch time of the CME in the corona. 

Similarly, for the HM model the arrival time calculation can be written as: 

, , , RiHM{A,di) 

taHM — toHM H 77 • (iDj 

VHM 

which reduces with Equation 13 to an expression consistent with MOll (their 
Equation B7). 

Note that for SSE, not only is the arrival time well defined, but the prediction 
if a spacecraft or planet is hit by an ICME or not is given by the conditions A > A 
(hit) and A < A (no hit). This is much more precise as compared to FP, where, 
strictly, a point never hits another point. Also, for HM, which corresponds to 
SSE for A = 90°, the circle will always hit if the in situ spacecraft is inside the 
half-space given by (/> ± 90°. 

2.2. Arrival speed 

Additionally, for an expanding circular front, every point along the front moves 
with a slower speed compared to the apex. Thus a prediction of a corrected 
arrival time should always include a corrected speed (see MOll). Because we 
assume radial propagation, the component of the velocity is always the radial 
component away from the Sun, and the others are zero. For self-similar ex- 
pansion, the following relation must be valid because the speed needs to be 
proportional to distance along the front, so the shape does not change with 
time: 

dj _ ViSSE ^^^^ 
RissE Vsse 
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Rearranging and substituting RissE from Equation 12 shows that the front at 
an angle A to the ICME apex moves with the foUowing speed {VissE < Vsse) 
away from the Sun: 



V.sse = Vsse'-^^^^^j^^^^^ (18) 
(1 + sm A) 

Again, for A = 90°, Equation 18 reduces to: 

V^HM = VhM COs(A) (19) 

as stated by MO 11 for the Harmonic Mean case. 



3. Plots for arrival time and speed corrections 

In this section we explore the significance of the previously derived formulas 
by plotting the corrections for arrival time and speed for spacecraft positions 
along the ICME flank, for different sets of parameters. In Figure 3 we introduce 
five SSE fronts with angular half widths of A = [10, 25, 40, 60, 90]°. These values 
correspond to a radius of curvature, which is equal to the radius r of the circle 
as given by Equation 4, of r = [0.086,0.211,0.321,0.433,0.5] AU for an apex 
distance of Rsse = 1 AU. We will also use Vsse = [400, 600, 800, 1000] km s-\ 
covering typical ICME speeds observed in situ at 1 AU throughout a solar cycle 
(Richardson and Cane. 2010). 

3.1. Arrival time 

In Figure 4a, wc plot the arrival time correction tcssE given by Equation 14 
as a function of A for an ICME speed of Vsse = 400 km s~^ and the in situ 
spacecraft at Earth distance {di = 1 AU). This is done similarly to MO 11 (their 
Figure 8b), but these authors restricted themselves to the HM case (A = 90°), 
which is the blue solid curve in Figure 4a. The functions for A = [10, 25, 40, 60]° 
each stop at the same value for A, though only seen in the plot for A=10°. 
This is because for A > A the ICME does not hit the spacecraft, and t^ssE 
becomes imaginary. For smaller values of A, a few degrees difference in A can 
lead to large differences in the arrival time, whereas for large A values the curves 
becomes increasingly flatter. Also note that for A = 60° there is already almost 
no difference to HM, so the formulas mainly affect cases where small values of A 
are combined with a A close to A, which is - in other words - the case for glancing 
encounters of ICMEs which are of small angular extent along the line-of-sight. 
This is either the case for intrinsically narrow ICMEs or elongated flux-rope like 
ICMEs which have a high inclination of their axis to the ecliptic plane. 

Figures 4b, 4c, 4d show similar plots for Vsse = [600,800,1000] km s^\ 
which allow to quickly look up, also useful for real-time assessments, how long 
the flank might be delayed for a given set of ICME parameters. As an example, 
it can be seen that for an average width of an ICME shock front (a full width 
of 80° corresponding to A = 40°), the arrival time correction stays low (< 10 
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hours) for encounters within A < 25°, but it quickly rises when A > 25° and 
may easily delay the front arrival by one or even two days. Thus, the conclusion 
is that the arrival time for an ICME front approximated by SSE is increasingly 
hard to pinpoint at its flanks, where small errors in A may lead to large errors 
in the arrival time prediction. 




10 20 30 40 50 o 10 20 30 40 50 



A [degrees] A [degrees] 

Figure 4. (a); The arrival time correction t^sSE assuming the ICME parameters 
VssE = 400 km s~^and the in situ observer situated at = 1 AU, for different values of the 
ICME half width A. (b) Similar to above, but for VssE = 600 km s"\ (c) VssE = 800 km s~\ 
(d) VssE = 1000 km s"!. 



3.2. Arrival speed 

For Figure 5 we repeated this analysis for the speed corrections. We use the 
same sets of parameters as in the previous section to plug into Equation 18. Note 
that Viss E does not depend on dt , because the front always moves with constant 
speed, regardless of the heliospheric distance of /. Again, where A > A the ICME 
does not hit the spacecraft. Clearly, the functional forms are governed by the 
same correction factor as compared to the arrival time in the previous section. 
Again, these plots are meant to be useful for looking up the effect of a glancing 
encounter for given parameters. For cases where A > 30°, the arrival speed can 
be slower in the order of 100 — 200 km s~^ compared to the apex. Again, this 
shows that for geometrically consistent calculations of arrival times and speeds it 
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is necessary to include these formulas because they lead to significant deviations 
from those inferred from the ICME apex. 




Figure 5. (a): The corrected ICME speed VissE a^t a longitudinal separation A between the 
ICME apex and the in situ target, assuming the ICME parameters VssE = 400 km s~^and 
the given angular half widths A. (b) Similar to above, but for VssE = 600 km s~^, (c) for 
VssE = 800 km s-l, (d) for VssE = 1000 km s^i. 



4. Conclusions 



Our aim was to show how arrival times and speeds can be analytically calculated 
for self-similar expanding circular fronts with constant angular width. We have 
shown plots of the corrections for speeds and arrival times for an observer at 1 AU 
and typical ICME speeds at 1 AU. This is not only useful for techniques which 
assume this geometry for application to observations by a single-spacecraft he- 
liospheric imager (HI) instrument (SSEF, Davies et al. 2012), or two His (Lugaz 
et al., 2010), to forecast the effects of interplanetary coronal mass ejections at 
Earth and other planets, but for any empirical/analytical technique which is used 
to calculate arrival times of ICMEs (summary by Siscoe and Schwenn, 2006; see 
also Vrsnak et al. 2012, this issue, for the drag-based model). Assuming the SSE 
geometry leads to the following results: For glancing encounters of an ICME 
with a planet or spacecraft, the flank parts of an ICME might be delayed by one 
or - in extreme cases - two days, and speeds can be affected in the order of 100 
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to 200 km s^^. These corrections are both significant, given that arrival time 
predictions are usually in the order of ±12 hours (e.g. Davis et al. 2011), and 
geoeff'ects of ICMEs arc thought to be determined by the reconnection electric 
field —Vr X depending on the plasma radial velocity component Vj. and the 
southward magnetic field component B-^ (e.g. Burton, McPherron, and Russell, 
1975). 

Further assumptions we used for the derived corrections arc that after the 
front passed the heliocentric distance of the planet or spacecraft, (1) the ICME 
propagates radially away from the Sun and (2) it does so with constant speed. 
Clearly, in a structured solar wind, the ICME front might easily become distorted 
by the heliospheric current sheet, co-rotating interaction regions or other ICMEs, 
thus one always has to keep in mind that the geometry used here is a highly 
idealized one. We also note that in the SSE framework, the curvature of the front 
is coupled to the full width, leading to a rather rigid description of its shape. De- 
coupling the width from the curvature would be possible in an elliptical model 
(or a related geometry used by Savani et al. 2011), at the cost of introducing 
another free parameter. We leave these ideas for future studies. 

Nevertheless, the expressions presented in this paper allow to check the results 
of SSEF against ICME in situ observations. It is indeed our aim to use the 
presented formulas in connection with the SSEF technique and its limiting cases 
(EPF and HMF) to test the efficacy of using heliospheric imagers as a predic- 
tion tool for real time space weather forecasting by studying a larger number 
of events. To do so in a geometrically consistent way, the derived corrections 
will be applied to comparisons between ICME arrival times and speeds derived 
from STEREO/HI with SSEF and those observed in situ by various spacecraft, 
extending other studies using elongation-fitting techniques and their stereoscopic 
versions (Davis et al. 2009; Liu et al. 2010a,b; 2011; Mostl et al. 2009, 2010, 2011; 
Lugaz et al. 2012). 
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